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We study the properties of strongly interacting Bose gases at the density and temperature regime 
when the three-body recombination rate is substantially reduced. In this regime, one can have a 
Bose gas with all particles in scattering states (i.e. the "upper branch") with little loss even at 
unitarity over the duration of the experiment. We show that because of bosonic enhancement, pair 
formation is shifted to the atomic side of the original resonance (where scattering length a s < 0), 
opposite to the fermionic case. In a trap, a repulsive Bose gas remains mechanically stable when 
brought across resonance to the atomic side until it reaches a critical scattering length a* < 0. 
For a s < aj, the density consists of a core of upper branch bosons surrounded by an outer layer 
of equilibrium branch. The conditions of low three-body recombination require that the particle 
number N < 1.024(T/oj) 5,/2 in a harmonic trap with frequency tj. 



One of the most fascinating aspects of quantum gases 
is the role of strong interaction, which is generated by 
bringing the system close to a Feshbach resonancepQ. 
The scattering between particles near resonance leads 
to a very large negative or positive s-wave scattering 
length a s , causing strong attraction or repulsion between 
atoms in scattering states. In the case of fermions, 
preparing the system on the atomic side of the resonance 
(a s < 0), the ground state of the system exhibits BCS- 
BEC crossover as the system is brought across the res- 
onance to the molecular side (a s > 0)j2j. In contrast, a 
repulsive Fermi gas is metastable, as fermions with pos- 
itive scattering length can form bound states (Feshbach 
molecules) through three-body recombination. Once Fes- 
hbach molecules are formed, they can decay into deep 
bound states through collisions, which leads to atom loss. 

In the two few years, there have been active experimen- 
tal studies on strongly repulsive fermions, driven by the 
quest of Stoner ferromagnetism[3]. While it is now shown 
that ferromagnetism is absent in repulsive fermions [1], 
many earlier experiments performed over a wide range 
of physical conditions all show similar behavior in atom 
loss (5][6j. Moreover, an early experiment [7] has clearly 
demonstrated that the energy density exhibits a maxi- 
mum near resonance on the molecular side of the reso- 
nance, which is found to be caused by Pauli blocking on 
the formation of bound pairs JB . 

In the case of Bose gases, attractive interactions will 
cause mechanical instability at low temperatures. Thus, 
most studies focus on repulsive Bose gases. However, 
like repulsive Fermi gases, repulsive Bose gases are only 
metastable. For weak interactions, the collision rate 
due to three-body recombination (73 = — n dn/dt) 
is 73 = c(47r/ia s /m)n(r£a 3 ), where c is a dimensionless 
constant [S], while the two-body collision rate is 72 = 
na 2 s v where v is the typical velocity of the bosons. For 
weak repulsion, n l ' 3 a a <C 1, 73 is sufficiently low that 
the system is essentially free of molecules. In the last 
two years, there are increasing number of experiments 
on strongly repulsive Bose gases at low temperatures (TUJ- 



il2j . However, at low temperatures, 73 increases rapidly 
in the strongly repulsive regime, i.e. n x ^a s > 1. This 
leads to severe atom loss as the system approaches res- 
onance, and the system is far from equilibrium. While 
one can explore strong interaction effects by bringing the 
system quickly in and out of the strongly interacting re- 
gion, it is not clear how to define equilibrium properties 
in such situations. 

The situation is different at higher temperatures and 
lower densities, i.e. lower fugacities. At temperatures 
T > T c , we have v ~ ySfesT/m ~ h/(m\), where 
T c is the the BEC transition temperature, and A = 
h/y/2TtmkBT is the thermal wavelength. Close to uni- 
tarity, a s in 72 and 73 is replaced by A in this tem- 
perature regime, and we have 72 = (fcsT/fi)(nA 3 ), and 
73 = C{k B T/h)(nX 3 ) 2 , where C = 9y/3/n ~ 4.96Q3]. As 
density drops, 73 will eventually fall below 72. And in 
the presence of a trap, the spatially averaged rate of total 
particle loss, — N~ 1 dN/dt = J dr7 3 n/ J dm = (73) a «ej 
will fall below the trap frequency uj, where (..} a ve means 
spatial average. 

In the density and temperature regime (referred to as 
"low-recombination" regime) where 73 < 72, (73)aue < 
lj, or 

nA 3 « 1, nA 3 < C~ 1/2 \/ 'tuu /k B T ', (1) 

where n 2 = J n? / J n, very few molecules are formed even 
at unitarity during the time when the Bose gas reaches 
global equilibrium through two-body collisions. We can 
then reach an equilibrium state where the bosons are 
in scattering states even though the system can accom- 
modate Feshbach molecules. This "low-recombination" 
regime has recently been realized by Salomon's group at 
ENS[ll]. In this paper, we shall point out a number of 
surprising properties of strongly interacting Bose gases 
in this low-recombination regime. We find that 
(I) Bose statistics enhances pair formation. As a result, 
molecule formation in a homogenous Bose gas is shifted 
to the atomic side (a s < 0), in contrast to fermions where 
the shift is to the molecular side (a s > 0) due to Pauli 
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blocking 8 . The energy change of the system when mak- 
ing transition from the upper to lower branch (defined 
later) is substantial even at temperatures as high as 10T C . 

(II) In a trap, when a repulsive Bose gas is brought 
across resonance in the low-recombination regime, it re- 
mains stable even on the atomic side (a s < 0), but up 
to a critical value a* < 0, and its density consists of a 
metastable "upper branch" core surrounded by an outer 
layer of bosons in thermodynamic equilibrium. Both re- 
gions are molecule free. The system will suffer mechani- 
cal instability for a s > a*. 

(III) The conditions for low-recombination at unitar- 
ity (Eq.Q) and mechanical stability constrain the total 
number of particles in a trap. In order to observe the phe- 
nomena in (II), we need N < aN* , N* = (k B T/hwf/ 2 
where the constant a = 1.024. For an estimate at 
T = 1/iK, lu = 2 7 r(250)sec~ 1 , we have N* ~ 6.5 x 10 4 . 

(A) Homogenous upper branch Bose gas: We first 
study the homogenous repulsive Bose gases that are 
free of molecules. Such system will be referred to as 
"upper branch" Bose gas, and is a good approxima- 
tion of a Bose gas in the low-recombination regime. 
In contrast, the equilibrium state of a Bose gas con- 
sisting of both atoms and molecules will be referred 
to as the "lower branch" or "equilibrium branch". To 
study the upper branch Bose gas, we use a general- 
ized Nozieres Schmitt-Rinkf NSR) [TB"] method recently 
developed by one of us (TLH) for the upper branch 
Fermi gas [8]. (We set both h and ^ to 1 from now 
on). It is straightforward to see that the equation of 
state is identical to that of a Fermi gas, except that 
all the Fermi functions are replaced by the Bose dis- 
tribution functions 71b(w) = l/(e"/ T — 1). The result 
is n(n,T) = n a ((i,T) + An sc (T,/j,) + An 6d (T,/z), where 
n (T,p) = X)k n s(£k) 1S the density of the ideal Bose 
gas; An sc (T,[i) and An bd (T,p) are the interaction con- 
tributions of the scattering states and the bound states 
respectively, 

An (ftT) = --^ — n B (u)—- — , (2) 



A^(u,T) = -^n s M g ))°^, (3) 
q 

where to(q) = q 2 /4m — 2/i. £(q, uj) is the phase of the 
inverse T-matrix in a medium whose explicit expression 
is given in ref.[S]. It arises from the branch cut of the 
T-matrix (i.e. the scattering states). 0Jb(q) is the pole of 
the T-matrix (i.e. the bound states), and is the solution 
of the equation 

m 1 v-^ / 7(k;q) 1 \ 

k \ m m / 
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FIG. 1. The phase diagram of a homogeneous upper branch 
Bose gas with fixed density n: At the blue curve that sep- 
arates the upper and lower branch, the energy density un- 
dergoes a discontinuous jump as shown in Fig[2] The purple 
dashed curve represents a state with n = 0. In the "unstable" 
region, the number equation for chemical potential does not 
have a solution. 

where 7(k;q) = 1 + n B (^ q/2+k ) + n B (£ q/2 _ k ) describes 
the bosonic enhancement of the medium on pair fluctua- 
tions, £fc = £fc — fi, and eu — k 2 /2m. 

As pointed out in ref.[8], the NSR results Eqns.Q and 
^ reduce to the scattering state and bound state con- 
tributions in the rigorous virial expansion in the low fu- 
gacity limit. The equation of state for the upper branch 
at lower temperatures therefore corresponds to ignoring 
An bd , whereas that of the equilibrium state (the lower 
branch) includes both An sc and An bd ; 

nupper (T, m) = n (T, ft) + An sc (T, /z) (5) 
n eqm i(T, pi) = n a (T, «) + An sc {T, u) + An bd (T, «). (6) 

With this prescription, we can calculate all the thermo- 
dynamic properties of these branches [TB]. Our results are 
summarized in the next two sections. 

(B) Phase diagram of homogenous upper branch Bose 
gas: It is useful to define for bosons an analog of "Fermi" 
momentum and "Fermi" temperature as kp = (Q^n) 1 / 3 
and Tp = k F /2m; and Tp/T c — 2.3, where T c = 
3.3n 2 / 3 /m is the BEC transition temperature. The phase 
diagram of the upper branch Bose gas for fixed n is shown 
in Figure [l] The corresponding behavior of the energy 
density at T = 4Tp is shown in Figure [2] 

The dashed purple line in Fig. [l] corresponds to the 
state of zero compressibility k = 0, where k — dn/dp. 
While the decrease in n as a s > increases is similar to 
that of Fermi gas, this phase diagram differs from that of 
Fermi gas in a fundamental way, as the transition from 
upper to lower branch is shifted to the atomic side of the 
resonance, (a s < 0). See Fig [l] and Fig (2j In other words, 
the stability of the upper branch Bose gas extends into 
the atomic side. This shift is due to Bose statistics. To 
see this, we note that the condition for determining the 
emergence of a bound state in the medium (i.e. Eq.Q) 
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FIG. 2. The upper panel is the energy density across reso- 
nance at T = \Ef = 9.2T C : rescaled by energy E {T) of a 
noninteracting system at the same temperature. The jump 
represents a transition from the upper to the lower branch. 
Even at this high temperature, interaction energy and the 
jump are substantial fractions of the total energy. The lower 
panel shows a c (q) as a function of q for fixed n. The jump 
is due to the sudden change in [i as the system switches 
branches. 



can be rewritten as 
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"-g(^q/2+k) + "-B(Cq/2-k) 
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(7) 



(8) 

which is the condition for emergence of bound state in 
vacuum with an effective scattering length a e //(q, u). In 
vacuum, a e ff reduces to a s , and the bound state occurs 
when 



> 0. In a Bose medium, a bound pair with 
total momentum q occurs when a e ff(q, uj — i-u(q)) > 0, 



or when 
1 



< 



a c where 

4?T »i?(£q/2+k) + n B (£,q/2-k) , , 

( q ) ~ n ^ v <u - w 



In other words, if one approaches the resonance from 
the atomic side, a bound pair with total momentum q 
will emerge at scattering length a c (q) < 0, which is on 
the atomic side of the original resonance. The values 
of — l/a c (q) at different q's are shown in Fig{2] That 
— l/a c (q) reduces to as q — > oo is because the effect from 
the bosonic medium becomes less important for large q, 
as in fermion case[H]. The boundary between the equi- 
librium branch and the upper branch corresponds to the 
emergence of bound pairs with q = 0, and is given by 
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FIG. 3. Phase diagram in a trap at fixed temperature T 
and trap frequency to: /J, a represents the state of compress- 
ibility k — 0. fj,^ separates the equilibrium branch with 
the unstable region. The latter has k < 0. Each point on 
this diagram denotes a density profile of the Bose gas, with 
fj, being the chemical potential at the trap center. The den- 
sity profile can be generated by an upward vertical trajec- 
tory using LDA, (see text). The three horizontal lines with 
arrows are trajectories of a Bose gas across resonance into 
the atomic side at fixed fi. The termination point (denoted 
by a black dot) on each /i-trajectory indicates the critical 
scattering length a* for that fi. For a 3 < a, (a s > a,), 
the density profile is stable (unstable). Hence, the density 
profiles (a) with chemical is stable, whereas the densi- 
ties (b) and (c), with chemical potential ^3 are unstable. For 
T = IfiK, and lu = 2ty(250Hz), we find (-A/o s )* = 0.14, 0.18 
and 0.21 for the trajectories —([i/T)i = 3, —(n/T)2 = 2.5 and 
— (/j,/T)3 — 2 respectively. The corresponding particle num- 
bers are (iVi, N 2 , N 3 ) = (2.1, 3.4, 5.2) x 10 4 . All these systems 
satisfy the condition to be in the low-recombination regime, 
Eq.|l]) or equivalently Eq.(10l, as (^ave/^ = 0.14,0.35, and 
0.96 respectively, and JVi, N 2 , N 3 < N* , where N* = 6.5 x 10 4 
for this temperature and trap frequency. 



a^ 1 (q — 0, T, £t) = 0, where \i is constrained by the total 
density n. Finally, we note from Fig(2]that as one crosses 
the resonance from the molecular to the atomic side, the 
energy change at the boundary between upper and lower 
branch is substantial even at T = ATp = 9.2T C . 

(C) Upper branch Bose gas in a trap: In a trap, the 
density profile within local density approximation (LDA) 
is given by n(r) = n upper {iJL — V(r),T), where [i is the 
chemical potential at the center of the trap. A global 
view of the density profile can be obtained from the phase 
diagram in the (/x/T)-(A/a s ) plane, Figure [3j where A is 
the thermal wavelength. The density profile along a ra- 
dial direction starting from the trap center corresponds to 
a vertical line emerging from —/J,/T upward. A trapped 
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FIG. 4. The density profiles of an upper branch Bose gas in 
a trap : Fig. (a), (b), and (c) correspond to the densities (a), 
(b), and (c) in Fig. [3] Density (a) has an upper branch core 
and an outer shell of equilibrium state. The dashed purple 
curve is the density of the equilibrium branch, Eq.|6f , at the 
same T, a s , ui, and N. The densities (b) and (c) are unstable 
as they contain regions where An/A^i < 0. For density (a), the 
width of the unstable region becomes less than inter-particle 
spacing and is therefore non-existent. 



Bose gas is therefore specified by a point (— X/a s , —fj,/T) 
on this diagram. 

Figure [3] describes the behavior of the Bose gas as it 
is swept from the molecular to the atomic side. Three 
regions are found from the equation of state: (i) the 
equilibrium branch, (ii) the upper branch, and (iii) a re- 
gion of mechanical instability where k < 0. (n > for 
both (i) and (ii)). The boundary between (i) and (iii), 
and that between (ii) and (iii) will be denoted as /i^ (T) 
and n^ a \T) respectively. fi^ a '(T) is the boundary of zero 
compressibility, k — 0. fi^ b ' (T) is the boundary where 
bound pairs with zero momentum begin to form. 

The unstable region (iii) intervenes between branches 
(i) and (ii). Because of this intervention, any density 
profile whose center starting from the upper branch on 
the atomic side will pass through the unstable region, and 
is therefore unstable. The width of the unstable region is 
Ar = r a -r b , where r a and r b are given by fi a = fj,-V(r a ) 
and /i b = fi — V{ v b)- As seen in Fig[3| the difference fj, b — 
H a (and hence Ar) becomes very small close to resonance. 
Should Ar be less than inter-particle spacing, i.e. Ar < 
n (f)" 1 / 3 , which occurs at a critical ratio A/a* for given 
/i/T, r a and r b can be viewed as a single point f. In this 
case, the unstable region disappears. The critical ratio 
A/a* can be estimated by setting Ar = rt(r a ) _1//3 . (a* is 
a function of T and fj,). 

Thus, for an upper branch Bose gas characterized by 
the point (— A/a s , — /x/T) on this diagram, it will only be 
stable when — l/a s < — l/a*(T, /i), so that the unstable 



region with a width Ar in real space collapses to zero. 
The vertical line labelled (a) in Fig{3] represents such a 
density. Its density profile is shown in Fig (4^,, which con- 
sists of an upper branch inner core and an equilibrium 
branch outer layer, both of which are free of Feshbach 
molecules. Compared to the density profile of the lower 
branch (dashed line in Figj4^i which is close to Boltz- 
mann distribution) , one sees a discernible kink in the up- 
per branch density. When — A/a s exceeds — A/a*(T, /z), 
such as (6) and (c) in Fig, 3, the corresponding density 
profiles are unstable, for they will contain a region of 
negative compressibility as shown in Fig|4jD and[4j;. 

From the equation of states for both the upper branch 
and the equilibrium branch, Eq.([5]) and ([6]), one can de- 
termine the total number of particles once the chemical 
potential at the center is specified. We also find that total 
particle number N changes little with a s for given /j, and 
T. It is straightforward to show that N — J drn(r) has 
the general form N — A(T/uj) 3 , where A is a dimension- 
less number depending on (— fi/T, —\/a s ). We find that 
A < 1. This is expected, as the critical number of Bose- 
Einstein condensation in a harmonic trap with frequency 
lo is N bec = (O.^^iT/ufiU^ . On the other hand, for 
the trapped gas to be in the low-recombination regime, 
Eq.|l]) imposes constraints on the central density, and 
hence the total particle number N. To find an estimate 
of this constraint, we approximate the actual density (say 
that in Fig. 4a ) by the Boltzmann form, which then gives 
N ~ e M / T (T/w) 3 [T7]. Within the same approximation, 
we find the quantity n in Eq.|l]) to be n = 3 3 / 4 e M / T /A 3 , 
which then imply 

N < N* — a(T/uj) 2 - 5 , a = 3 3/4 C~ 1/2 = 1.024. (10) 

We have thus established the results (I) to (III) . 
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